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In three-dimensional turbulent flows energy is supplied at large scales and cascades down to the 
smallest scales where viscosity dominates. The flux of energy through scales implies the generation of 
small scales from larger ones, which is the fundamental reason for the irreversibility of the dynamics 
of turbulent flows. As we showed recently, this irreversibility manifests itself by an asymmetry of 
the probability distribution of the instantaneous power p of the forces acting on fluid elements. In 
particular, the third moment of p was found to be negative. Yet, a physical connection between the 
irreversibility manifested in the distribution of p and the energy flux or small-scale generation in 
turbulence has not been established. Here, with analytical calculations and support from numerical 
simulations of fully developed turbulence, we connect the asymmetry in the power distribution, i.e., 
the negative value of (p^), to the generation of small scales, or more precisely, to the amplification 
(stretching) of vorticity in turbulent flows. Our result is the first step towards a quantitative 
understanding of the origin of the irreversibility observed at the level of individual Lagrangian 
trajectories in turbulent flows. 


I. INTRODUCTION 

The generation of small scales, or large velocity gradi¬ 
ents, is one of the most striking physical phenomenon of 
3-dimensional (3D) turbulent fluid flows, and is respon¬ 
sible for a flux of energy e from large to small scales. 
Remarkably, in the limit of very small viscosity or very 
large Reynolds number, the third moment of the lon¬ 
gitudinal velocity difference between two points sepa¬ 
rated by a distance x, {Au(x)^), is related to the en¬ 
ergy flux by the relation (Au(x)^) = —^sx, which is one 
of the very few exact results in turbulence theory [l|. 
In elementary terms, two points are more likely to be 
pushed closer together (repelled) when their relative en¬ 
ergy is large (small) [^. This fundamental asymmetry 
persists all the way down to very small distances so the 
third moment of the velocity derivative d^Ux is negative: 
{{dxUx)^) < 0. In fact, available data from experiments 
using hot-wire anemometry Si and from direct numer¬ 
ical simulations (DNS) have led to the conclusion that 
the normalized third moment oi dxUx, i.e., the skewness, 
^ {{dxUxf)/{{dxUxff/\ is negative and approx¬ 
imately —0.5, with at most a weak dependence on the 
Reynolds number Si. In homogeneous isotropic flows, 
the seminal work of Betchov Q shows that the third mo¬ 
ment of dxUx is related to the generation of small scales in 
turbulence, through amplification of vorticity by vortex 
stretching. 

Because of the existence of an energy flux from large 
to small scales, turbulence is a non-equilibrium phe¬ 
nomenon, thus intrinsically irreversible. The possibil- 
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ity to probe turbulence by following the motion of indi¬ 
vidual particles in both numerical and laboratory high- 
Reynolds-number flows [I®: leads to new insights on 
irreversibility and offers new opportunities for quantita¬ 
tively understanding turbulence [uHil. 

Recently, we observed that the energy differences along 
particle trajectories present an intriguing asymmetry: ki¬ 
netic energy grows more slowly than it drops along a tra¬ 
jectory [^. The consequence of this asymmetry is that 
the third moment of the power p = a- u is negative, where 
u and a are the velocity and acceleration of the fluid (see 
[Tsil for a related discussion). As a possible explanation, 
one may expect the pressure gradient, which dominates 
the fluctuations of the power, to provide an explanation 
for the negative sign of the third moment of p (l6l - [l8l| . 
Unexpectedly, however, in 3D, the contribution of the 
pressure gradient to the third moment of power is very 
small (T^ . 

Here, we provide a physical relation between the nega¬ 
tive third moment of p and the generation of small scales 
by turbulence, i.e., vortex stretching. In the following, it 
is convenient to decompose the power as 

P = PL+PC ( 1 ) 

where pl = u-a^, = u-9tu and pc = u- ac = u- (u- V)u 
are the local and convective parts, respectively. We find 
that the magnitude oi p = pl A- pc is much smaller than 
the magnitudes of its components pl and pc, which im¬ 
plies significant cancellation between pc and pl- On av¬ 
erage, the magnitude of pc is larger than that of pc. Note 
however that the cancellation between p^ and pc does 
not automatically follow from the well-known cancella¬ 
tion between and ac since pl, Pc involve 

only one component of a^, ac. We demonstrate that 
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the moments of p, up to the third order, are dominated 
by the moments of pc- In particular, the third moment 
(p^) has the same sign as (pc)- We show analytically 
that (pq) is a surrogate for vorticity amplification. This, 
together with the observation that (pq) determines the 
sign of (p^), leads us to the conclusion that the origin of 
the negative sign of the third moment (p^) comes in fact 
from small scale generation, thus clearly establishing a 
relation between the generation of small scales and the 
observed irreversibility in the flow. 

II. NUMERICAL METHODS 

A. Direct Numerical Simulation of Navier-Stokes 
Turbulence 

We investigated numerically turbulent flows, obtained 
by solving directly the Navier-Stokes equations: 

5tu(x, t) + (u(x, t) ■ V)u(x, t) 

= —VP(x, t) + z/V^u(x, t) + f(x, t) (2) 
V-u(x,t)=0 (3) 

where u(x,t) denotes the Eulerian velocity field, P is 
the pressure, v is the viscosity, and f(x, t) is a forcing 
term; the mass density is arbitrarily set to unity. Solving 
the equations in a simple cubic box of size (27r)^ with 
periodic boundary conditions allows us to use efficient 
pseudo-spectral methods. 

The forcing term acts at large scales, or equivalently, 
on Fourier modes at low wavenumbers, |k| < Kf. It 
is adjusted according to a method proposed in [22|, in 
such a way that the injection rate of energy, e^, remains 
constant: 

fk = e, — -if |k| < Kf , (4) 

with Kf = 1.5. In the code units, the energy injection 
rate Ei has been set to Ei = 10“^. Note that in station¬ 
ary turbulent flows, the energy injection rate equals the 
energy dissipation rate, Ei = e. 

The code is fully dealiased, using the 2/3-rule 
method [2^. We have chosen two different resolutions, 
corresponding to the highest resolved wavenumber of 
^max = 256 and /cmax = 384 (effectively equivalent to 
768 and 1152 grid points in each spatial direction), with 
the corresponding values of the viscosity v = 1.6 x 10“^ 
and 9.0 x 10“®, respectively. With these values, the Kol¬ 
mogorov scale p = (y^ IeY^'^ is such that the product 
^max X p is very close to 2 in both cases, ensuring ad¬ 
equate spatial resolution. The corresponding Reynolds 
numbers are Rx ~ 193 and 275, respectively. 

Once expressed in terms of spatial modes, Eq. m re¬ 
duces to a large set of ordinary differential equations, 
which were integrated using the second-order Adams- 
Bashforth scheme. The time step St has been chosen so 
that the Courant number Co = Urms'kmaxSt ^0.1, where 


Rx 

193 

275 

430 

(/)/£= 

3.83 X 10^ 

7.36 X 10^ 

1.32 X 10® 

ipDle^ 

2.15 X 10^ 

5.00 X 10® 

1.20 X 10"* 

ipDle^ 

1.78 X 10® 

4.25 X 10® 

1.07 X lO'* 

-{pL Pc)/e^ 

1.77 X 10® 

4.26 X 10® 

1.07 X lO'* 

\^pl)/{e^Rl) 

0.87 

0.99 

0.96 

/3 

0.83 

0.86 

0.90 


TABLE I. Second moments of the distributions of p/e, 
pc/e and ph/e at the three Reynolds numbers studied in 
this article. The correlation coefficient between pc and pl 
is approaching —1 as Reynolds number increases. The val¬ 
ues of P are measured from fitting the conditional averages 
(Pl|pc> = -Ppc- 


Rx 

193 

275 

430 

-{P^)le^ 

3.87 X 10® 

1.23 X 10* 

3.21 X 10* 

-{Ph)/e^ 

5.39 X 10-* 

2.40 X 10® 

1.00 X 10® 

{PcPL)/e^ 

4.54 X 10-* 

2.05 X 10® 

8.99 X 10® 

-{pcpl)/e^ 

4.02 X 10'* 

1.84 X 10® 

8.29 X 10® 

{Pl)/e^ 

3.44 X 10* 

1.63 X 10® 

7.63 X 10® 

C = {pIp)/{pI) 

0.108 

0.088 

0.066 

l-P 

0.17 

0.14 

0.11 

{p^)/{ph) 

0.072 

0.051 

0.032 

1-/3-C 

0.061 

0.052 

0.037 

-{pI)/{pI) 

0.64 

0.68 

0.76 

P-2C 

0.62 

0.69 

0.77 

{pcpl)/{ph) 

0.75 

0.77 

0.83 

13-C 

0.72 

0.77 

0.83 

-{phPL)/{pl) 

0.84 

0.85 

0.90 


0.83 

0.86 

0.90 


TABLE 11. Third moments of the distributions of p/e, pc/e 
and pe/e at the three Reynolds numbers studied in this arti¬ 
cle. The last 8 rows compare the normalized moments with 
our predictions. 


Mrms is the root mean square value of one component of 
velocity. 


B. Data from the Johns Hopkins University 
Database 

We also used additional numerical simulation data at 
Rx = 430 from the Turbulence Database of the Johns 
Hopkins University. The flow is documented in [^. We 
computed the statistics presented here with at the mini¬ 
mum 2 X 10® points. 
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III. THEORETICAL BACKGROUND 
A. Elementary relations 

To investigate the moments of p, pc and we first 
note that pc reduces to a simple form that is particularly 
useful, namely: 

PC = u • (u • V)u = u • S • u, (5) 

where the rate of strain tensor S is the symmetric part 
of the velocity gradient tensor Vu: S = [Vu+ (Vu)^]/2 
or s., = {diUj + djUi)/2. Geometrically, the straining 
motion decomposes into a superposition of compression 
or stretching along three orthogonal directions, denoted 
by Bi, with three straining rates, A^. The vectors and 
the straining rates Aj are the eigenvectors and eigenval¬ 
ues of S. A positive (respectively negative) value of \i 
corresponds to stretching (respectively compression) in 
the direction e^. Volume conservation (incompressibil¬ 
ity) imposes that Ai -I- A 2 -f A 3 = 0, i.e., the amount of 
stretching and compression along the three directions 
sums up to 0. 

Equation ([5]) shows that in a steady (frozen) flow, the 
kinetic energy of a fluid element changes only through the 
action of the rate of strain. The antisymmetric part of 
the velocity gradient tensor expresses the local rotation 
rate in the fluid, and is characterized by the vorticity 
w. The expression for the amplification (stretching) of 
vorticity in the flow is given by (w • S • w) [ 2 ^ [ 2 ^. In 
a statistically homogeneous flowjthe following identity 
holds: (w • S • w) = —|(tr(S^)) Q. Last, we note that 
in a homogeneous isotropic flow, the second and third 
moments of dxUx can be simply expressed in terms of the 
moments of tr(S^): {{dxUx)^) = ^(tr(S^)), and tr(S^): 
{{.dxUxf) = i(tr(S3)) a. 


B. Decomposition of power p: order of magnitudes 


The magnitudes of the fluctuations of the convective 
and local components of p may be estimated from sim¬ 
ple dimensional arguments: \pc\ ~ \pl\ U'^/tk, where 
U is the typical size of the velocity fluctuations, and tk 
is the fastest time scales of the turbulent eddies. Using 
the known relation tk ~ {U^/e)/R\ [H, [ 2 ^, one finds 
\pc\ ~ \pl\ ~ where Rx is the Reynolds number 
based on the Taylor microscale, and characterizes the in¬ 
tensity of turbulence. The growth of the variances of 
Pc/s and Phis, as R \, as predicted by this simple dimen¬ 
sional argument, is found to be consistent with our DNS 
results, see Table HI This result sharply contrasts with the 
fact that the variance of p is known to grow more slowly 
with the Reynolds number, as R^^^ 141. This difference 
in the observed scalings as a function of the Reynolds 
number is due to a very strong cancellation between pp 
and Pc, see Table HI We observe that the magnitudes of 
the third moments (pcPp), with m -\- n = 3, are found 


to increase with m, see Table HD signaling that the con¬ 
tribution of Pc to the third moments is more significant 
than that of pp. In fact, as we will show, the sign of (p^) 
is dominated by (pc)- 

Although the cancellation between pc and pc is remi¬ 
niscent of the well-documented cancellation between ac 
and ap [312 jwe stress that it cannot be deduced from 
the results of [3 13 ■ In fact, pc and pp involve the 
projections along the direction of the velocity u, of ac 
and ap, respectively. Our results therefore show that the 
cancellation between ac and ap affects their components 
along the velocity direction, which does not result auto¬ 
matically from [3112 ■ 

In the following subsections, we begin by expressing 
(pc) in terms of vortex stretching, before establishing 
the prevalence of pc on (p^). 

IV. RESULTS 

A. Vortex stretching and moments of pc 

It is convenient to express pc (Eq. ([5|)) by project¬ 
ing the velocity u and the rate of strain S in the basis 
of the three perpendicular unit vectors characterizing 
the straining motion. In this basis, the velocity u is de¬ 
composed as: u = UiGi, where = u • is the 
coordinate of u along the direction e^, and the rate of 
strain tensor is expressed as S = AiB^Bi. Denot¬ 

ing Xi the cosines of the angles between the velocity u 
and the unit vectors Bp Xi = u • Bi/|u| = Mi/|u|, the 
expression of pc reduces to: 

PC = ^ = 

i=l i=l 

In a turbulent velocity field, small wave numbers (or 
large scales) provide the main contribution to the veloc¬ 
ity held, u, whereas the rate of strain S is determined by 
the large wave numbers (or small scales). The two helds 
u and S are therefore expected to be only weakly corre¬ 
lated. Let us now assume that S and u are uncorrelated. 
This approximation implies that the three cosines, Xi, are 
uniformly distributed between —I and I. Geometrically, 
the three cosines are the coordinates of a point that is 
uniformly distributed on the unit sphere in 3D. This as¬ 
sumption allows us to compute the averages necessary to 
evaluate explicitly the third moment of pc- 

Namely, Eq. ([5]) leads to: 



= y (7) 

The assumption that u and S are uncorrelated also im¬ 
plies that all the cosines Xi, {i = I... 3), are inde¬ 
pendent of the eigenvalues of S. As a consequence. 
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(A™a;") = (A™)(i") for any m and n. Using the observa¬ 
tion that {xi,X2,X3) represents the coordinates of a point 
that is uniformly distributed on the unit sphere, which 
gives the symmetry relations such as (x®) = (x®) = (x®), 
(xfx|) = {x^xf) = (x|xi), etc, we therefore obtain 

+i{x\xl)l ^ A^aA 

-f 6 (xiX 2 X 3 ) (A 1 A 2 A 3 ). ( 8 ) 

The averages of the products of Xi in Eq. (I5|) can be calcu¬ 
lated by using elementary geometrical considerations Q 
and the results are: 

{xf) = y , {x\xl) = ^ , (x^x^x^) = As' ( 9 ) 

Substituting Eqs. ([HI) and dH) into Eq. (ED, and using the 
incompressibility of the flow, Ai -I- A 2 -I- A 3 = 0, leads to 
the following expression for the third moment of pc' 

{pD = A(|u|6)(AiA2A3). (10) 

Using the relation (A 1 A 2 A 3 ) = (l/3)(tr(S^)) = —(l/4)(a;- 
s-u;) dig, one finally obtains: 

{pD = A(|u|6)(tr(S3)) = -|(|ur)(u; • S • cu). (11) 

Thus, Equation (El) relates the third moment of pc 
to vortex stretching in such a way that positive vortex 
stretching ((w • S • w) >0) gives rise to a negative value 
of {pI). 

Furthermore, many experimental and numerical stud¬ 
ies show that the probability distributions of individual 
components of velocity u are close to Gaussian, with 
small deviations that can be quantitatively explained (see 
e.g., [H, [ 2 ^). Assuming a Gaussian distribution of u al¬ 
lows us to express the 6 *^’' moment of velocity in Eq. ED 
in terms of the velocity variance (u^). Using other known 
identities, in particular concerning the relation between 
(tr(S^)) and the skewness of the velocity derivative ; 
as explained in Appendix A, Eq. ED can be written as: 

(pI) = ^So^uxRle^- ( 12 ) 

The weak dependence of the velocity derivative skewness 
on the Reynolds number, Sg^u^ oc ii in Eq. ED 
suggests a small correction to the simple order of mag¬ 
nitude analysis for the third moment: {pc) oc with 
(5 ft: 0.1. 

The assumptions of lack of correlation between u and 
S, and of a Gaussian distribution of the velocity u, also 
lead to an exact determination of the variance of pc- 
(Pc) ~ Appendix A. This expression for the 


second moments of pc provides further justification for 
the dimensional estimate of the variance of pc, and is 
found to be in very good agreement with our DNS results 
(see Table [ID. 

Having established the relation between the third mo¬ 
ment (pc) ^ad vortex stretching, we now establish that 
the third moment (p^) is dominated by (pc)- To this 
end, we first consider the cancellation between pc and 
PL- 



FIG. 1. The joint probability density function (PDF) between 
pc/s (horizontal) and px/e (vertical) at R\ = 275, color- 
coded in a logarithmic scale (see color-bar). Equal-probability 
contours, separated by factors of 10, are shown. The PDF is 
concentrated close to the pc + Pl = 0 line, indicating that 
the two quantities pc and px are nearly anti-correlated with 
each other. The black dashed line shows {px\pc)/ e, which is 
approximately —0.86 x pels. The white dashed line shows 
{pc\px)/ e, which is approximately —pxjE. 


B. Cancellation between px and pc 

We note that in homogeneous and stationary flows, the 
first moments of p, pc and px are all exactly 0. Table |T| 
shows that the correlation coefficient between px and pc'- 
{pcPl)I{{Pc){pI)V^'^, is approximately —0.9 and seems 
to approach —1 as the Reynolds number increases. This 
strong anti-correlation results in significant cancellation 
between pc and px, so the variance of p is much smaller 
than those of pc and px- 

Although the range of values of Rx covered by the 
present study is not sufficient to reach unambiguous con¬ 
clusions, our results are generally consistent with the ex¬ 
pected scalings: (p\) ~ ip^) oc R\, and (p^) oc [T^. 

Further insight into the strong cancellation between pc 
and Px can be gained by studying the joint probability 
density function (PDF) of pc and px, shown in Fig.|T]for 
our flow at Rx = 275. Fig. [1] clearly indicates that with a 
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FIG. 2. (a) The joint PDF between pc/e (horizontal) and p/e (vertical), and (b) The joint PDF between pn/e (horizontal) 

and p/e (vertical), all at R\ — 275, color-coded in a logarithmic scale (see color-bar). The equal-probability contours shown 
are separated by factors of 10. Note that the ranges of values of p/e shown are smaller than those of pc/s and plIe. In (a), 
the black dashed line represents the conditional average of p/e on pc js, and is very close to (p|pc)/e ~ 0.14 x pc/ e- The white 
dashed line corresponds to the averaged of pc /e, conditioned on p/e, and is well approximated by (pc'|p)/e « p/e. In (b), the 
black dashed line represents the conditional average of p/e onpe/e, and is very close to (p|p_L)/e « 0. The white dashed line 
indicates to the average of pe/e, conditioned on p/e, and is also very well approximated by (pl |p)/£ ~ 0. 


high probability, the values of pc and pp are concentrated 
close to the line pc +Pl = 0, thus implying a significant 
cancellation between the two quantities. The observed 
tendency of the correlation coefficient between pc and 
Pl to approach —1 as the Reynolds number increases 
implies that the joint PDF otpc andpp becomes increas¬ 
ingly concentrated around the line pp+pc = 0 al higher 
Reynolds numbers. In all our numerical simulations, we 
find an approximately linear relation between the con¬ 
ditional average (pplpc) and pc (shown as the black 
dashed line in Fig. [T]): {pl\pc) ~ —P{R\)pc, where the 
dimensionless coefficient (3{Rx) depends weakly on Rx- 
In agreement with the observed tendency of pp and pc 
to become increasingly anti-correlated as Rx increases, 
we find that P{Rx) slightly increases with the Reynolds 
number, see Tabled] This implies that {p\pc) ~ 
where the coefficient 1 —/3 decreases as Rx increases, from 
Ri 0.17 at Rx = 193 to ~ 0.10 at Rx = 430. We also ob¬ 
serve that the average of pc conditioned on pp, shown as 
the white dashed line in Fig. [1] is almost exactly equal to 
—Pi, which implies that (p|pl) « 0. 


Fig. Hjshows the joint PDFs of pc and p (a) and of pp 
and p (b). The conditional averages (p|pc) and (p|pl) 
are shown as black dashed lines, whereas the conditional 
averages {pc\p) and (pl|p) are shown as white dashed 
lines. The conditional averages of pc and pp on p have 
the particularly simple forms: (pc|p) ~ P and (pi|p) ~ 
0. In addition, the joint PDF of p and pp is almost 
symmetrical to both p = 0 and pp = 0. The power p is 
therefore well correlated with pc, but not with pp. 


C. Prevalence of pc on the moments of p 

1. General assumptions 

The prevalence of pc on the statistical properties of p 
shown by our numerical results leads to the conclusion 
that the second and third moment of p are expressible 
in terms of the corresponding moments of pc- To justify 
this claim, we use the two following results. 

A The numerical results shown in Fig. [1] demonstrate 
that at a fixed value of pc, Pl\pc fluctuates around the 
mean value (pp|pc) ~ —Ppc- This immediately implies 
the following relations: 

(PLPC) = -/3(Tc) and (ppp^) = -PiPc)^ (13) 

which can be easily justified by writing pp conditioned 
on a value of pc as: 

Pl\pc = -Ppc + i\pc, (14) 

where ^|pc is a random variable with zero mean and its 
distribution depends on pc- Eq. (fl^ is found to be nu¬ 
merically extremely well satisfied, see Table HH as a direct 
consequence of the quality of the linear dependence be¬ 
tween (ppIpc) and pc- 

B The lack of correlation between p on pp, demon¬ 
strated in Fig. |2| and manifested by the two relations 
(pp|p) ~ 0 and (p|pp) ^ 0, implies that: 

(ppp) Ri (p^pp) Ri (ppi) Ri 0. (15) 

The equalities shown in Eq. (HU) are only approximate. 
In the following, we explore the consequences of the in¬ 
dependence between p and pp by assuming for now that 
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these equalities are exactly satisfied, leaving for later a 
discussion of the errors made. 

As we show below the approximations A and B above 
lead to a very accurate prediction of all the second mo¬ 
ments of Pc, Pl and p, in terms of (p^) and /3. The 
predictions concerning the third moments, however, are 
not as accurate as those concerning the second moments, 
as a consequence of quantitative deviations from the sym¬ 
metry assumption B. 

2. Second moments 

Using Eq. (US]), ipcPc) = -l3{Pc)^ and Eq. (HH), 
{pPl) = 0 , we determine the second moment of pc as 
a function of (p^): (p|,) = l3{pc), from which we obtain: 

(p2) = (l-/?)(p^). (16) 

We find that the condition of independence (ppl) = 0 
is very well satisfied, which implies that Eq. (USD is nu¬ 
merically very accurately satisfied (see Appendix B and 
Table HVl). 

3. Third moments 

The results from Eq. (fT^ . showing that (p^) ~ —e^R\, 
together with the observation that (p^) ft! —e^R\ |14l| . 
also point to a strong cancellation between pc and pc 
in the third moment (p^). To relate the properties of 
the third moments of p to {p%), we begin by noting that 
Eq. dTSl) leads to the following expressions for the third 
order moments: {p\) = —I3{p^) and (plPc) = l^{p%), 
and hence to the expression (p^) = (1 — I3){pq). These 
expressions predict simple relations between the various 
moments (pcPl.) with m -|- n = 3 and (p^), and lead to 
the correct sign of (p^), thus justifying our claim that the 
assumption of independence beween p and pc imposes 
that the sign of (p^) is given by (pq). 

The expressions obtained above, however, are quanti¬ 
tatively not accurate. The reason is that while {plP^) is 
found to be very small (of the order of 1 % of |(Pc)l)i the 
numerical values of {p\p) are found to be much larger, 
of the order of 10% of |(Pc)|. The small, but significant 
error in {p\p) = 0 therefore leads to a significant reduc¬ 
tion of the numerical value of (p^), consistent with the 
numerical values shown in Table IE To take the effect of 
non-zero (p|p) into account, we denote C, = (pip)/(Pc)i 
where C is a positive number of order ^0.1 (see Ta¬ 
ble IJ). and decreases when R\ increases. This then leads 
to {pipe) = {P- 0{Pc) and (pi) = -{13 - 2C,){p%), and 
consequently: 

(p3) = (l-/3-C)(pi:). (17) 

Using Eq. (HZl) and the relation between (p^.) and vor¬ 
tex stretching, Eq. dnD, we obtain: 

(/) = -^(1 - /? - C)(|u|®)(w • S • w), (18) 


which establishes a quantitative relation between the 
time irreversibility, as measured by (p^), and vortex 
stretching, a small-scale generation mechanism in 3D tur¬ 
bulence. 

We note that Eq. (lEI), together with the observed 
scaling (p^) oc R^^^ and (p^) oc suggests that 

(1 — /?) oc Similarly, the dependence (p^) oc i?i, 

together with the observation of that (p^) ft —e^R\, 
imply, using Eq. (fT7)l . that (1 — /3 — C) oc R'^^ or 
1 — t^/{l — (3) oc ■ These are consistent with the 

values obtained numerically. As shown in Table IE the 
value of 1 — Q/{1 — (3) decreases slightly, from 0.36 to 0.34, 
when the Reynolds number increases from R\ = 193 to 
430. 

The results presented here thus show that, while the 
statement of independence between p and pi is merely 
an enticing approximation, taking quantitatively into ac¬ 
count the deviations from Eq. cni) does not affect our 
main conclusion: the third moment of p is controlled by 
the third moment of pc- 



FIG. 3. The PDFs of the cosines of the angles between the 
direction of the velocity u, and the eigenvectors ei of S at 
Rx = 275. The cosine of the angles is given by the inner 
product of the two unit vectors: Xi = Si-eu, where Bu = u/|u|. 
A lack of correlation between u and S would lead to a flat PDF 
of Gi ■ e„: U(|ei ■ e„|) = 1, which are supported by the data 
to a good degree. The small departure from this expectation 
points to a weak correlation between u and S. 


D. Lack of correlation between u and S 

We return here briefly to discuss the essential assump¬ 
tion that u and S are uncorrelated. Specifically, we exam¬ 
ine in this subsection the correlation between the angles 
of u and the eigenvectors and between the magnitude 
of u and the eigenvalues Xi. In the following, the values 
Xi are sorted in decreasing order: Ai > A 2 > A 3 . 
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FIG. 4. The thick lines show the conditional averages of 
the eigenvalues of S, (Ai|u^), made dimensionless by tk = 
1/(2(S^))^/^ = The thin lines with the same 

color are the corresponding unconditional averages (Ai), all 
at R\ = 275. For small values of u^, the conditional aver¬ 
ages are nearly the same as the corresponding unconditional 
averages, consistent with the assumption that u and S are 
uncorrelated. For large u^, the magnitudes of the conditional 
averages increase with u^, indicating a deviation from the as¬ 
sumption. Note that because the probability of having large 
values of decreases very rapidly with u^, the increases of 
(Ai|u^) at large have only very small effects on the uncon¬ 
ditional averages Ai. 


Fig. [3] shows the PDFs of |xi| = je^ • e„|, the absolute 
value of the cosine of the angle between the eigenvec¬ 
tor ei and the unit vector in the direction of the veloc¬ 
ity e„ = u/|u| (the sign of this cosine is immaterial) at 
R\ = 275. A complete lack of correlation between u and 
S implies that the PDFs of |ei-e„| are constant and equal 
to 1. Fig. [3] shows that this is close to be true. Namely, 
the probability of alignment between ei and u, i.e., of 
|e„ • ei| being close to 1, is slightly reduced. On the con¬ 
trary, the probability of alignment between e„ and 62 is 
slightly increased. The deviations observed numerically 
are weak, less than ~ 10 %, compared to the uniform dis¬ 
tribution. The cosine between and 63 , is very close 
to being uniformly distributed. The nearly uniform PDF 
of |ei • 011 1 indicate that the assumption that S and u 
are uncorrelated, explicitly used in the determination of 
(Pc)’ provides a very good first order approximation. 

The assumption that u and S are uncorrelated also 
implies that the conditional averages of the properties of 
S should be independent of the magnitude of u. 

Fig. |4] shows that the dependence of the conditional 
average of the eigenvalues of S on u^, (Ai|u^), is weak. 
Systematic deviations are visible at large values of u^, 
where the magnitudes of the averaged conditional eigen¬ 
values are larger. The probability of large values of u^, 
however, drops very rapidly when increases (28l - [30l| . so 
the effect of this weak dependence of Ai on has only 
a small effect on the low-order moments of pc studied 


here. 

In summary, the results presented here and in the Ap¬ 
pendix C show that the assumption of a lack of corre¬ 
lation between u and S provides a very good first-order 
approximation to describe the third moment of pc- 


V. DISCUSSION AND CONCLUSION 

Our work, aimed at understanding the third moment 
of the power p acting on fluid particles (p^) « —£^R\ , 
and its implication for the physics of turbulent flows [l4| , 
rests on decomposing p into two parts: a local part, 
Pi, = u • 9tU, induced by the change of the kinetic en¬ 
ergy at a fixed spatial point, and a convective part, 
PC = \i- V(u^/2), due to the change in kinetic energy 
along particle trajectories, assuming the velocity field is 
frozen. We observe that the two terms pc and pi cancel 
each other to a large extent, resulting in a much smaller 
variance of p compared to those of either pc or pi. This 
cancellation may be qualitatively explained by invoking 
a fast sweeping of the small scales of the flow by the large 
scales [ 1 ^. In physical words, kinetic energy along par¬ 
ticle trajectories, is mostly carried (swept) by the flow, 
and changes far less than it would change by keeping the 
flow fixed, or by varying the flow with the same position 
in time. This fact has been documented in a slightly dif¬ 
ferent context [Slj. Our results provide a quantitative 
characterization of how much sweeping reduce the indi¬ 
vidual contributions of pi and pc- 

One of the two main results of our work is that the 
third moment of pc, expressed in terms of the rate of 
strain, S, and the velocity, u, as: pc = u • S • u, can be 
exactly determined, by using the physically justified ap¬ 
proximation that u and S are uncorrelated. Remarkably, 
we find that (p^) is directly related to vortex stretching, 
(w • S • w). In particular, the negative sign of (p^) origi¬ 
nates from the positive sign of the vortex stretching, due 
to small-scale generation by turbulence. This observa¬ 
tion provides the first basis for our claim that the third 
moment of p is related to the generation of small scales 
in 3D turbulent flows. 

The other main observation of our work is that, despite 
the strong cancellation between pi and pc, the power 
p correlates with pc, but not with pi, as revealed by 
the nearly vanishing conditional averages (p|pi) ~ 0 and 
(pi|p) ~ 0. Assuming these conditional averages are ex¬ 
actly zero leads to a simple relation between (p^) and 
(Pc)- The (weak) corrections to this simple assumption 
modify only quantitatively the results. 

Taken together, these two observations, namely that 
(p^) is controlled by (p^) and that (p^) is directly linked 
to vortex stretching, (w • S • w) >0, allow us to establish 
a relation between the third moment of power, (p^), and 
vortex stretching. Thus, the recently observed manifes¬ 
tation of irreversibility in studying the statistics of indi¬ 
vidual Lagrangian trajectories can be understood as re¬ 
sulting from small-scale generation in 3D turbulent flow. 













For lack of essential information concerning the quanti¬ 
ties investigated here, our work rests on several assump¬ 
tions supported by numerical observations. The well- 
known fact that velocity, u, and the rate of strain, S, are 
dominated by large- and small-scales, respectively, makes 
it plausible that these two quantities are mostly uncor¬ 
related. Our numerical results confirm this expectation. 
Although small, and of little relevance for the low-order 
moments studied here, the deviations observed suggest 
an interesting structure, which would be worth elucidat¬ 
ing. The observation that p and pl are not correlated, in 
the sense that the conditional averages {p\pl) and (pl|p) 
are both very close to zero, rests only on numerical ob¬ 
servations, and requires a proper explanation. Under¬ 
standing and quantifying the weakness of the correlation 
between p and pl in 3D turbulent flows may provide im¬ 
portant hints not only on higher moments of p, but more 
importantly, on the structure of the flow itself. 

We note that studying the cancellation between pc 
and PL by directly focusing on the effect of the pressure 
gradient, — u • VP is likely to lead to satisfactory results 
when studying the second moments of p, as the pressure 
term has been documented to provide the largest con¬ 
tribution to the variance (p^) In 3D, however, the 

third moment — ((u- VP)^) has been shown to contribute 
negligibly to (p^), whose understanding requires the in¬ 
vestigation of other correlations [l^ . 

Finally, the arguments provided here to explain the 
negative third moments of power fluctuations of particles 
in 3D turbulent flows should not be applied to 2D turbu¬ 
lence, in which the third moment of p is also negative, and 
grows with a similar power of the Reynolds number 1^ , 
but the amplification of large velocity gradients is due to 
entirely different physical processes (32|. Still, one may 
expect that the manifestations of irreversibility, in 2D 
turbulence as well as in a broad class of non-equilibrium 
systems, to be fundamentally related to a flux in the sys¬ 
tem. 
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APPENDIX A: ALTERNATIVE EXPRESSIONS 
OF THE THIRD MOMENT OF pc 

In many experiments and numerical simulations, the 
probability distributions of individual components of ve¬ 
locity u are close to Gaussian [2^. That observation 
allows us to estimate explicitly the 6*^ moment of |u| in 


Rx 

193 

275 

430 

225{pl)/i7e^Rl) 

VTS y q 

~T~ ^ ^vc 

-0.24 

-0.30 

-0.37 

-0.38 

-0.40 

-0.42 


-0.52 

-0.62 

-0.67 


TABLE III. Third moments of the distributions of pc/e and 
p/e at the three Reynolds numbers studied in this article. 


terms of the second moment: 

(|u|4) = ^(|u|Y, (|ur) = |(|u|y, (19) 

which, when substituted into Eq. (ED, gives an expres¬ 
sion for {p^) as: 

(pl) = ^(|up)3(tr(S3)). (20) 

The same assumptions and similar elementary alge¬ 
braic manipulations as those used to establish Eq. (ED, 
lead to an exact expression for the second moment of pc '■ 

{Pc) = ^(|u|^)(tr(S2)) = (21) 

This relation is found to be in very good agreement with 
our numerical results, see Table HI 

The known relation between tr(S^) and the experimen¬ 
tally accessible moments of dxUx 0- 

10 '! 

= ( 22 ) 

together with the expression for the second moment of 

dxUx' {{dxUx)"^) = e/(15z/) [2^,[2^, leads to the following 
expression for the third moment of pc' 

(Pc) = (23) 

where = {{dxUxY)/{{dxUx)'^)^('^ is the skewness of 

the velocity derivative. 

Combining Eq. (121) and (l23l) gives the following rela¬ 
tion between the skewness of dxUx and the skewness of 
PC- 


Sa^u. = (24) 

The values of as well as the ratio 

225{p^)/{7e^R\), determined from our numerical 
simulations, are shown in Table uni The correponding 
value of the skewness of dxUx, Sa^u^ is found to be 
approximately Ri —0.4, which is well within the range 
of values of Sa^u^ reported from experiments and 
simulations at comparable Reynolds numbers 
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R\ 

193 

275 

430 

{PL Pc)/{pi)R^{pl)^^^ 

-0.907 

-0.923 

-0.944 

_/31/2 

-0.912 

-0.928 

-0.947 

{p^)/{ph) 

0.178 

0.147 

0.110 

1-d 

0.17 

0.14 

0.10 

(pi)/(pc) 

0.823 

0.851 

0.892 

/3 

0.83 

0.86 

0.90 


TABLE IV. Parametrisation of the second moments of pc/£, 
Pl/s and p/e, compared to the expression in terms of P given 
by Eq. I25I26I and 1271 at the three Reynolds numbers studied 
in this article. 


APPENDIX B: PREVALENCE OF pc ON THE 
MOMENTS OF p 

A. Second moments of p 

The decomposition of the distribution of pL condi¬ 
tioned on Pc, see Eq. da. involving a random variable 
P/\pc with a zero mean, together with the assumption of 
independence betweenp andpi, (ppi) = 0 (see Eq. (Ilbp l. 
leads to a full description of the second moments of pc, 
Pi in terms of j3 only: 

(pi) =/3(Pc) (25) 

{PLPC)= -P{pI) = (26) 

(p2) =il-l3){pl) (27) 

The values of P are measured from the conditional aver¬ 
age (plIpc) = —Ppc, Eq. (fT4)l - This allows us to check 
how accurately are Eqs. (l25l) to (l27ll satisfied. The nu¬ 
merical results are shown in Table HVl The relations given 
by Eqs. (ESI) to dSil) are found to be very well satisfied, 
with very small errors, thus demonstrating that the pro¬ 
posed parametrization in terms of {pq) and /3 provides a 
very good description of the second moments of p. 


B. Conditional averages {p\pc) and (pc|p) 

Crucial to the argument relating the third moment (p^) 
to the third moment of pc is the observation that the con¬ 
ditional averages (p|pl) and (pl|p) are close to 0. Fig¬ 
ure 2b of the main text shows these averages (pl, {p\pl)) 
and ((pl|p),p), which appear as horizontal and vertical 
straight lines respectively on the scale of the figure. Our 
argument is then based on identities such as: 

/ OO 

'P{pL){p\PL)pldpL, (28) 

-OO 

where P{pl) is the PDF of pL. Equation (l28)l shows that 
if {pl\p) = (pIpl) = 0, then, (p^pi) = (pip) = 0. 

Possible deviations from zero of the moments (pj^p) 
and {p^pl) therefore indicate that the conditional aver¬ 


ages (pIpl) and {pl\p) are not exactly zero. These mo¬ 
ments can be readily estimated from the various third 
moments {pcp'l) with m-|-n = 3 shown in Table 2 of the 
main text: 

(p^Pi) = {PcPl) + 2(pcpi) + (pi) (29) 

and 

(PlP) = (pi) + (pipe). (30) 

When compared to (pq), the value of {p^pi) is approx¬ 
imately zero ( |(p^Pl)/(pc)I < 2 %), but {pip)/{pl/) < 
10%. This points to a departure of (p|pl) from being 0, 
which we explore here. 

Figure [5] shows the conditional averages of (p|pl) for 
the three direct numerical simulation (DNS) runs dis¬ 
cussed in this article, and also the integrand in Eq. (I28p . 
For all three cases, the curves differ weakly but consis¬ 
tently from 0. While for p^ < 0, the values of (p|pl) are 
very small, they differ noticeably from 0 on the positive 
Pl side, with an approximately linear dependence on p^. 
In order to examine the effect of this deviation of (p|pl) 
on (^ = {pip)/{pc), we non-dimensionalize the variables 
PL and p in Fig.[5]by (pp)^^^- In particular, in Fig. [51(b), 
we plot 'P{pL){p\PL)pi/{Pc)^^'^ ■ In this way, the areas 
under the curves in Fig.jSKb) give (pip)/(pc)^'^^ = C^pey 
where Sp^ = {Pc)/{Pc)^^^ is the skewness of pc, which 
depends weakly on the Reynolds number as shown in Ta¬ 
ble 2 of the main text. Fig.jHb) shows that (p/p) / {pc)^^"^ 
decreases when the Reynolds number increases. This is 
consistent with the observed decrease of the value of C, 
with the Reynolds number. 

As shown in Table 2 of the main text, / = 0.11, 0.088 
and 0.066 at R\ = 193, 275, and 430, respectively. In 
fact, over the limited range of Reynolds number that we 
studied here, we observed that C/(l — P) remains well 
below unity: //{I — jd) ~ 0.64, which ensures that the 
third moment of p, as given by Eq. [12] in the main text, 
is determined by (p^): 

(p3) = (l-/3-C)(p^) = (1-/3) (l - {pl). (31) 

We note that the scaling of (p^) ~ R\ reported be¬ 
fore M, together with the scalings {pD ~ R^ and 
(1 — /3) ~ obtained in this work, implies that 

1 ~ [C/(l ~ /)] ~ These predictions can only 

be checked by using DNS at much higher Reynolds num¬ 
bers, and with adequate statistical resolution. 

For comparison. Fig. | 6 | shows the conditional average 
of Pi on p: (pijp) and the integrand of (p^Pi). As it was 
the case in Fig. |5l the two quantities p and pl are nor¬ 
malized by {pc)^^^ in the way such that the areas under 
the curves in Fig. | 6 l(b) represent the normalized moment 
{p^Pl)I{pI)'^I'^ = {{p^PL)/{Pc))Spc- A systematic devi¬ 
ation of (pijp) from being 0 is visible in Fig. EKa). On 
the other hand, the integrand 'P{p)p^{pl\p) is noticeably 
non-zero only in a small range of p, which results in a 
much smaller value of {p^pi) compared to {pip). 
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FIG. 5. The conditional average oip on pl- (a) raw data {p\pl) / vs. PlUpc)^^'^ and (b) 'P{pl){p\pl)Pl/ (Pc)^^^, which 
provides the integrand in Eq. (I28II to calculate the normalized third moment = C'S'pc- conditional average 

{p\pl) differs weakly, but consistently from being 0, especially for pL > 0, where the conditional average (p|pn) decreases 
approximately linearly. This leads to an appreciable negative contribution to (pip): C = (pip)/(Pc) ^ 10% in the range of 
Reynolds number studied. 




FIG. 6. The conditional PDF of pL on p: (a) (pi:,|p)/(pc)^^^ vs. pUpcV^^, and (b) 'P{p){pl\p)p^ /■ The areas under the 
curves in panel (b) give the normalized third moment (p^pn)/(pc)^^^i which are much smaller than those for (pip)/(pc)^^^ in 
Fig.IHb). 


APPENDIX C: LACK OF CORRELATION 
BETWEEN u AND S 


The results presented in the main text give a good in¬ 
dication that assuming u and S are uncorrelated provides 
an appropriate first-order approximation. This is in fact 
corroborated by the quantitative agreement between the 
numerical results and the predictions. Here we present 
further information concerning the correlations between 
u and S. 


The first hint of a correlation between the velocity and 
strain was provided by Fig. 3 of the main text, which 
showed weak, but noticeable deviations from a uniform 
distribution for PDFs of the cosines of the angles be¬ 
tween the direction of velocity, e^, and the eigenvectors 
of strain, especially between e„ and ei and between 
and 02 . 

It may be expected that the alignment between u and 
the eigenvectors of S depends on the magnitude of u. 
Figure [7] shows the joint probability distribution func¬ 
tion between and je^ • e^j for i = 1 (a), f = 2 (b) and 
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FIG. 7. Dependence of the alignment between e; and on at Rx = 275. Panels (a) (respectively (b) and (c)) show 
P(u^, |eu • Gil), the joint PDF of and |e„ ■ ei| for i = 1 (respectively i = 2 and i = 3). The indicated colour coding refers to 
the decimal logarithm of the PDF, i.e., log 7^(0^, |e„ • ei|). The equal-probability contours, shown as dashed lines, are close 
to, but deviate from being vertical, which indicates a systematic dependence on u^, especially for i = 1 and f = 2. Panel (d) 
shows the conditional average (|ei ■ e„||u^) vs. (also shown as dash-dotted lines in Panels a-c), which weakly deviates from 
being constant as implied by the assumption of lack of correlation between S and u. 


1 = 3 (c). The bending of the equal-probability contours, 
shown as dashed lines in Fig. [7Ka-c) reveals a weak, but 
systematic dependence of the alignment between u and 
as a function of u^. Fig. [71(c) shows that the effect 
is considerably weaker for z = 3 than it is for i = 1 and 
i = 2. We observe that the average of |ei-e„| conditioned 
on u^, plotted as the dash-dotted lines in Fig. jT]) a-c) and 
separately in Fig. IlKd), shows clear variations as a func¬ 
tion of u^, especially for z = 1 and z = 2. The informa¬ 
tion presented in Fig. [7] thus reveals that influences 
not only the eigenvalues of strain, as shown in Figure 4 of 
the main text, but also the statistical properties concern¬ 
ing the orientation of u with respect to the eigenvectors 
of S. This in turn induces a dependence of the third and 
fifth moments of pc on that is more complicated than 


the expectation based on the lack of correlation between 

u and S. 

The results shown by Fig. [7] thus show small, but sys¬ 
tematic deviations of the statistical quantities relevant to 
Pc from the expected dependence on . In comparison, 
the dependence on the magnitude of S seems to be much 
weaker. Figure [8] shows that the value of the average of 
the cosines of the angles, |e„ • e^j, conditioned on |S|, de¬ 
pends significantly less on |S|: the variations shown in 
Fig. [8] are of the order of ~ 5%, whereas the ones shown 
in Fig. |71i are of the order of ^ 20%. This difference 
points to a stronger dependence of the alignment prop¬ 
erties on the large scale features of the flow, than on the 
small scales. 

The results presented in this section thus demonstrate 
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that, while the results obtained in this work by assuming 
that u and S are uncorrelated do provide a good approx¬ 
imation to the third moment of pc, small, but system¬ 
atic deviations from this assumption are visible. Judging 
from the present results, the dependence on |u| seems to 
be generally more important than the dependence on | S |. 


FIG. 8. Dependence of the alignment between ei and e„ on 
|S| for i = 1 (full line), i = 2 (dashed line) and i = 3 (dash- 
dotted line) at Rx = 275. The variation of the conditional 
average of |ei -eul on |S| is much weaker than the dependence 
on shown in Fig. [T] 
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